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ABSTRACT 

Strong lensing is to provide accurate mass measurements within Einstein radii. However, precise modelling of strong lensing systeins 
can be affected by external mass distributions, e.g. the group or cluster within which the lens is embedded. In this article, we propose to 
turn this limitation to our advantage and to use precise strong lensing modelling to probe external mass distributions surrounding the 
lens. We consider SL2S J08544-0121, a galaxy group at z = 0.35 that contains a strong lensing system. A simple elliptical isothermal 
potential cannot reproduce satisfactorily the strong lensing constraints. We include an external mass perturbation corresponding to 
the group within which the lens is embedded. The lensing properties of this perturbation are parametrised by its total mass M and a 
smoothing scale ^ that quantifies the characteristic scale over which M is distributed. For a range of these parameters, we are able to 
reproduce accurately the observations. This suggests that light is a good tracer of mass. Interestingly, this also shows that a localised 
strong lensing analysis (on scales of ~ 10") allows us to constrain global properties of the group as a whole (on scales of ~ 100"). 
Indeed, we constrain the group mass-to-light ratio to be M/L = 98 ± 27 (/ band, solar units, not corrected for evolution) and .v = 20"± 9 
(2 cr confidence level). We demonstrate that these strong lensing only constraints are due to the perturbed strong lensing configuration, 
where the main arc is located at ~ 5" from the galaxy, whereas its counter-image is found at ~ 8". To test independently our resulting 
strong lensing model, we pursue an independent weak lensing analysis of the group and find a mass-to-light ratio in the range 66-146 
(1-cr confidence level). We conclude that strong lensing can be used to probe mass distributions beyond their Einstein radius, and we 
characterise the kind of perturbed strong lensing systems allowing such an analysis: a non dominant strong lensing system used as 
a particle test probing the main potential. This kind of analysis needs to be validated further on other systems and may constitute a 
quick method to probe mass distributions. This is particularly relevant in the context of forthcoming wide field surveys which will 
yield thousands of strong lenses, some of them being perturbed enough to pursue the analysis proposed in this paper. 
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produced at TERAPIX and the Canadian Astronomy Data Centre as 
part of the Canada-France-Hawaii Telescope Legacy Survey, a collab- 
orative project of NRC and CNRS. Also based on HST data, program 
10876 and Keck telescope data. 



1. Introduction 

Gravitational lensing allows one to measure line-of-sight pro- 
jected mass distributions. When the surface mass density of a 
lens is larger than a critical threshold, the light from a back- 
ground source galaxy is lensed into multiple images, a regime 
called strong lensing (SL). These multiple images provide sttong 
observational constraints on the projected mass disttibution of 
the lens within the Einstein radius. Since the discovery of the 
first gravitational arc in the core of massive galaxy clusters 
twenty years ago (Lynds & Petrosian 1986; Soucail et al. 1987), 
strong lensing is widely used to probe mass disttibution of sttuc- 
tures on different scales: galaxies (see, e.g. the SLACS sur- 
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vey, Koopmans et al. 2006), galaxy clusters (see, e.g. Halkola 
et al. 2006) and recently galaxy groups (Cabanac et al. 2007; 
Limousin et al. 2009; Belokurov et al. 2009). 

1.1. Perturbing a Strong Lensing System 

Because most of the galaxies in the Universe are part of larger 
structures, either groups or clusters, many SL systems are also 
part of larger structures, (see, e.g. Kundic et al. 1997; Fassnacht 
& Lubin 2002; Faure et al. 2004; Morgan et al. 2005; Williams 
et al. 2006; Momcheva et al. 2006; Auger et al. 2007; Tu et al. 
2008; Auger et al. 2008; Grillo et al. 2008; Treu et al. 2009; 
Inada et al. 2009). Any mass at a small angular distance of a 
strong lens may induce measurable perturbations in the lens- 
ing signal. Not considering this external perturbation can seri- 
ously bias the results inferred from the SL modelling as shown 
by Keeton & Zabludoff (2004): they found that SL modelling of 
dual-image lenses, which neglect environment, grossly overes- 
timate both the ellipticity of the lens galaxy (Ae/e ~ 0.5) and 
the Hubble constant (Ah/h ~ 0.22). Besides, standard models of 
foiu'-image lenses, which approximate the environment as a tidal 
shear, recover the ellipticity reasonably well but still overesti- 
mate the Hubble constant (Ah/h ~ 0.15). They argue that most 
of the biases are caused by neglecting convergence from nearby 
massive neighbours. More generally, the topic of modelling a 
lens with an extemal mass perturbation has been addressed by 
different authors (see, e.g Keeton et al. 1997; Kochanek et al. 
2001; Keeton & Zabludoff 2004; Oguri et al. 2005, 2008; Dye 
et al. 2007, and references therein). 

To summarise, a precise SL modelling can be affected by 
some extemal mass distribution, and people have tried to take 
this bias into account in order to improve the SL modelling. 
As observations become more and more accurate, we can ex- 
pect to be more and more sensitive to external mass distributions 
near strong lenses. In this article, we propose to turn this bias 
to our advantage and to use the perturbations measured in SL 
modelling as probes of the external mass distribution. First, we 
remind the reader of such an attempt to probe a galaxy cluster 
potential using a perturbed SL system. 

1.2. The "Ring" Test in Abell 1689 

In the core of galaxy cluster Abell 1689, Limousin et al. (2007) 
reported SL systems ("rings") formed around three elliptical 
galaxies located 100" away from the cluster centre, i.e. between 
the strong and weak lensing regime regions. These SL systems 
should be sensitive to the extemal shear and convergence pro- 
duced by their parent cluster (Kochanek & Blandford 1991). 
Based on simulations, Tu et al. (2008) showed that such strong 
lenses could be used to probe the cluster potential locally. Tu 
et al. (2008) applied this method to the three rings discovered 
in Abell 1689. They found that modelling these three rings only 
(i.e. without including any other multiply-imaged systems also 
produced by the cluster) provides strong evidence for bimodality 
of the cluster core; it is not possible to model simultaneously the 
three rings assuming a unimodal mass distribution for the clus- 
ter This bimodality confirms previous parametric SL studies of 
Abell 1689 (Miralda-Escude & Babul 1995; Halkola et al. 2006; 
Limousin et al. 2007; Leonard et al. 2007; Saha et al. 2007; 
Okura et al. 2008). More importantly, this result shows that SL 
features of l-2"-wide Einstein rings actually contain informa- 
tion on the mass distribution of the parent cluster, i.e. on a much 
larger scale than their Einstein radii. In other words, this study 



suggests that strong lenses can be used to probe mass distribu- 
tions beyond their Einstein radius. In this article, we develop this 
idea on another perturbed SL system located in a galaxy group, 
SL2SJ08544-0121. 

All results are scaled to the flat, ACDM cosmology with 
£1m=0.3, Qa=0.7 and a Hubble constant Ho=70 kms ' Mpc~'. 
In such cosmology, at z = 0.35, 1" corresponds to physical trans- 
verse distance 4.94 kpc . All images are aligned with the WCS 
coordinates, i.e. north is up, east is left. Magnitudes are given in 
the AB system. Ellipticities are expressed as (a^ - b^)l{a^ + b^), 
and position angles are given counter-clock wise with respect to 
the west. Luminosities are given for the i band, in solar units, 
not corrected for passive evolution. Shear and convergence are 
computed for a source redshift z,=1.268. 

2. SL2SJ08544-01 21: Presentation & Data 

SL2S J08544-0121 is part of the Strong Lensing Legacy Survey 
(SL2S, Cabanac et al. 2007) , which collects SL systems in the 
Canada France Hawaii Telescope Legacy Survey (CFHTLS). 
SL2S J08544-0121 is a galaxy group at z = 0.35 presented by 
Limousin et al. (2009) which contains a SL system (Fig. 1). 

2.1. Ground-Based Imaging 

SL2S J08544-0121 has been observed in five bands as part of 
the CFHTLS. The j'-band data are used to build luminosity maps 
from isophotal magnitudes of elliptical group members and to 
perform a weak-lensing analysis. 

Fig. 1 bottom panel shows a 10' x 10' CFHTLS i-band im- 
age, The white cross gives the location of the strong lens. We 
draw luminosity isodensity contours of 10^, 3x10^, lO*", 3x10^ 
and 10^ L0kpc"^. Fig. 1 top-right panel also shows a CFHTLS 
l-arcmin^ gri colour image centred on the lens. 

2.2. Space-Based Imaging 

The strong lensing features detected from ground-based images 
have been observed with the Hubble Space Telescope (HST). 
Observations were done in snapshot mode (CIS, P.l. Kneib, ID 
10876) in three bands with the ACS camera (F814, F606 and 
F475). Fig. 1 shows a colour image of the strong lens based 
on these observations. We report two multiply-imaged systems: 
the first system is bright and forms a typical cusp configuration 
perturbed by a satellite galaxy (labelled Dwarf on Fig. 1). The 
second system is a very faint arc located West of the lens at a 
larger radius. It is not possible to reliably separate or identify in- 
dividual images on the faint arc that may be singly imaged only. 
Moreover, given its faintness, spectroscopy is hopeless with cur- 
rent facilities, as the surface brightness is ca. 31 magarcsec"^. 
Therefore, this faint arc is not used in the following analysis. 
As can be appreciated on Fig. 1, the HST data brings significant 
additional information on the lensed features. 

2.3. Spectroscopy 

We have used the Low Resolution Imager and Spectrograph 
(LRIS, Oke et al. 1995) on the Keck telescope to measure the 
spectroscopic redshift of both the lens and the brightest arc of 
the SL2SJ08544-0121 system. On January 14 2007, we ob- 
tained 300 seconds on the lensing galaxy and 4 exposures of 900 
seconds on the arc, using a 1.0" width slit. A 600 lines mm"' 
grism blazed at 4000 A and a 400 lines mm ' grating blazed at 
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8500 A were used in the blue and red channels of the instru- 
ment, both light paths being separated by a dichroic at 5600 A. 
The corresponding dispersions are 0.6/1.85 A and resolutions 
are 4.0/6.5 A in the blue/red channel. The resulting extracted 
spectra are shown Fig. 2. The lens presents a typical elliptical 
spectrum at z = 0.3530 ± 0.0005 with prominent H and K Call 
lines, 4000 A break, and G band. The spectrum of the arc shows 
a strong emission line at 8454 A, resolved in a doublet separated 
by ~ 5 A . This gives us an unambiguous redshift determination 
ofz = 1.2680 ±0.0003. 



3. Modelling the Lens 

In this Section, we attempt to reproduce the SL multiple images 
using a single elliptical isothermal potential centred on the bright 
galaxy. All optimisations are performed in the image plane us- 
ing the Lenstool software (Julio et al. 2007). We quantify the 
goodness of the fit by using the image plane RMS and the cor- 
responding x^- When necessary, we compare the fits using the 
Bayesian Evidence. 



2.4. Global Properties 



3.1. Observational Constraints 



The luminosity contours of SL2S J08544-0121 are elongated in 
the east-west direction. The SL deflector is populated by a single 
bright galaxy. Note (Fig. 1) that the innermost luminosity iso- 
density contour at 10^ L0 kpc"^ encompasses the SL system but 
also two bright galaxies located ~ 54" east from the SL system, 
making this light distribution bimodal. This is the only group of 
the sample presented by Limousin et al. (2009) for which the 
luminosity isodensity distribution is not clearly dominated by 
the lens, making this configuration rather exceptional: the large 
Einstein radius (~ 5") points towards a significant mass concen- 
tration associated with this lens, but the luminosity isodensity 
distribution is actually bimodal. 
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Fig. 2. 1-D spectra of the lens (up) and of the bright arc (low). The green 
shaded region masks the residuals of a strong OH atmospheric emission 



As explained in Section 2.2, we do not use tiie faint arc in the 
analysis and focus on the bright multiply-imaged system. This 
system is composed of 4 main images: 1.1, 1.2, 1.3 and 1.4. 
Additional images are produced by the satellite (dwarf) galaxy. 
These are not considered in the analysis because we do not want 
to complicate the modelling by adding a sub-halo for the satellite 
galaxy. We also make the assumption that, given its small size, 
the satellite galaxy does not influence images 1.1 to 1.4. That 
assumption will be discussed in Section 8 and addressed further 
in a forthcoming publication. Since the merging arc composed 
by images 1.2 and 1.3 is well resolved, we can safely associate 
two other images on this arc, namely 2.2 and 2.3 (Fig. 1). Their 
counter-images expected near images 1 . 1 and 1 .4 are not safely 
identified, therefore we do not use them in the analysis. Indeed, 
parametric strong lensing analyses are highly sensitive to images 
misidentifications and we prefer to use only images we are con- 
fident in. This gives us a total of 8 observational constraints. 



3.2. Shape of the Bright Galaxy 

In this subsection, we describe the properties of the light dis- 
tribution of the bright galaxy populating the strong lensing de- 
flector. We use the IRAF task ellipse to measure the shape of its 
isophotes. We find an ellipticity of e = 0.206 and a position angle 
of 39+5 degrees at 2" from the centre. A closer inspection of the 
galaxy centre clearly reveals a double core even though the outer 
isophotes are elliptical. The above measurements therefore cor- 
respond to the superposition of the light from each component. 
The spectrum of the galaxy presented in Fig. 2 do not show fea- 
tures of another galaxy at a different lens redshift along the line 
of sight. Given the similar colours of these two components, the 
bright galaxy may be the result of a recent merger. 



3.3. Positional Uncertainty 

The sizes of the multiple images have been estimated using the 
IRAF task imexamine. They range from 0.11" to 0.15", with 
a mean of 0.13". Therefore, the positional uncertainty is set to 
0.13". We note that this positional uncertainty sounds large for 
HST data but we stress that we are not in the case of point- 
like objects (quasar lensing) where the astrometric precision can 
reach 0.01". In our case, the images are extended and the depth 
of the snapshot observations does not allow to resolve better the 
conjugated points. We note that large positional uncertainties are 
often used in the case of extended images (see, e.g. Oguri 2010). 
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Fig. 1. Group SL2S J08544-0121 at z,,pec = 0.353. Upper Left: composite CFHTLS gri colour image (1 arcmin^ = 297 x 297 kpc"). Upper Right: 
composite HST/ACS F814W-F606W-F475W colour image (24" x 24" = 118 x 118 kpc^). We show the proposed multiple-image identification. 
The dwarf galaxy and the main extra image it produces are labelled. Lower: CFHTLS ; band (10' x 10' = 2969 x 2 969 kpc"). Luminosity 
isodensity contours of 10', 3 10', 10*, 3 10* and 10' L0 kpc - are drawn (continuous black line), and the white cross shows the location of the SL 
system. 
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3.4. Mass Model 

The lens potential is parametrised by a dual Pseudo Isothermal 
ElUptical Mass Distribution (dPIE, see Ehasdottir et al. 2007). 
The 3D density distribution of the dPIE is: 

Po 

P(r)=- , , 2 TTTT; '■cut>'-core- (1) 

(1 + r2/4re)(l + r^/rjj 

This distribution represents a spherical system with scale radius 
rcut, core radius rcore and central density po. 

This profile is formally the same as the Pseudo Isothermal 
Mass Distribution (PIEMD) described in Limousin et al. (2005). 
Its scale radius is set to 250 kpc, i.e. larger than the range where 
the observational constraints are found. Allowing rco,e to vary 
produces models with core radii much smaller than the radius 
over which we have observational constraints. Therefore, we can 
set rcore=0, and the dPIE profile becomes close to isothermal 
in the constrained region. The remaining free parameters of the 
dPIE profile are: 

- the halo centre position (X,Y) allowed to vary within 3" of 
the light distribution centre 

- the halo ellipticity e, forced to be smaller than 0.6, as sug- 
gested by numerical simulations (Jing & Suto 2002) 

- its position angle allowed to vary between and 1 80 de- 
grees 

- The fiducial velocity dispersion' allowed to vary between 
200 and 900 km/s. 

We emphasise that this fiducial velocity dispersion is not the 
Spherical Isothermal Sphere velocity dispersion. It is usually 
smaller, see EHasdottir et al. (2007) for a self-contained descrip- 
tion of the dPIE profile. 

3.5. Results: Bad Optimisation 

Results of the optimisation are given in Table 1. The fit is bad, 
the RMS error of image positions is 0.38" in the image plane 
(i.e. significantly larger than the assumed positional uncertainty 
of 0.13") and a reduced of 29. The halo position is found 
coincident to the light distribution centre within error bars. The 
halo ellipticity is at the higher edge of the input prior, and the 
position angle is equal to ~ 18 deg. Only when we allow the 
halo ellipticity to reach values as high as 0.9, are we able to 
reproduce the observational constraints (RMS equals to 0.06" 
for e = 0.9, 6~ 19.5 degrees, and the halo centre is offset from 
the light distribution centre by one arcsecond. 

We conclude that a single potential does not satisfactorily re- 
produce the observational constraints. In the rest of the paper, we 
add the contribution of the external mass distribution in the lens 
modelling. Finally, we have used the non-parametric method of 
Suyu et al. (2006) in parallel to our method and found that the 
observational constraints used in this work require an external 
shear component in order to be properly reproduced. 

4. An External Mass Perturbation Based on the 
Light Distribution: Does Light traces IVIass? 

The large scale properties of SL2SJ08544-0 121 shown on Fig. 1 
together with the failed modelling attempted in the last Section 
suggest the need to take into account an external mass perturba- 
tion. In order to test the hypothesis that light traces mass, this 

' linked to the central density by: cr^ = iGnpn -, 3 



external perturbation will be mapped from the known light dis- 
tribution properties. 

4.1. Luminosity Maps 

The first step is to build luminosity maps of SL2S J08544-0121 
from which we will derive the external mass perturbation prop- 
erties. Group members have been colour selected. All galax- 
ies having a r - i colour difference smaller than 0.15 with the 
bright galaxy deflector are selected as group members (Limousin 
et al. 2009). Because we want to describe the perturbation of the 
galaxy group on the SL system, this luminosity map should not 
take into account the light coming from the galaxies populating 
the SL deflector. Therefore, we select all group members except- 
ing the bright galaxy populating the deflector and the associated 
satellite galaxy. This partial group luminosity is referred to as 
Lext in the following. From this catalogue, we generate luminos- 
ity maps as described below. An important ingredient of this pro- 
cedure is the smoothing scale. Indeed, the shape of the luminos- 
ity maps strongly depends on the smoothing scale. Mechanically 
that has an influence on the properties of the derived mass maps. 
Therefore, this smoothing scale is a free parameter describing 
the external mass perturbation. 

We use the following smoothing scheme: the 10' x 10' 
CFHTLS /-band image is divided into cells of size c pixels, 
which translates into cx 0.186 ". We compute the rest- frame i- 
band luminosity L of each galaxy located in a given cell with 

L = io(*^G"-*^+'^'^+''^/^-^ (2) 

where M is the /-band isophotal magnitude of the galaxy, Mq 
is the solar absolute magnitude in the ; band; DM is the dis- 
tance modulus, and li the k-correction factor, estimated from 
elliptical templates by Bruzual & Chariot (2003) using single- 
burst stellar formation models. Then we sum up the luminosi- 
ties of all galaxies per cell to get the total luminosity within 
each cell. The resulting luminosity isodensity map is then con- 
volved with a Gaussian kernel of sigma w. This gives an angu- 
lar smoothing scale s equals to cx 0.186 xw" (since the pixel 
size equals 0.186"). Fig. 3 shows three luminosity maps, where 
we distribute the same total luminosity Lext for three smooth- 
ing scales. We draw luminosity isodensity contours of 10^ and 
10^ Lq arcsec"^. One can appreciate how the smoothing scale s 
influences the shape of the luminosity isodensity contours. 

4.2. From Light Map to Mass Map 

Once a luminosity map with a given smoothing scale s is ob- 
tained, we assume a mass-to-hght ratio Mext/Lext = 1 to convert 
it into a mass map. Because we omit the bright galaxy populating 
the deflector when building the luminosity map, this mass map 
can be considered as the external mass perturbation. Therefore 
we refer to the mass contained in this map as Mext- Then, we use 
the algorithm developed by Julio & Kneib (2009) to interpolate 
this mass map into a grid of analytic circular dPIE potentials, 
supported by Lenstool. 

We model flie mass distribution of SL2S J08544-0121 witii 
a 5' hexagonal grid of dPIE potentials. In order to build an 
adaptive grid where resolution follows the 2-D mass density, 
we recursively split the input mass map into equilateral trian- 
gles until the mean surface density per triangle is lower than 10^ 
M0 arcsec"^. Then we place a dPIE potential at each node of the 
grid with the following parameters: core radii rcore equal to the 
local grid resolution and cut-off radii rcut = 3 x rcore- In Julio 
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(To (km s ') 


RMS 


x' 


log(Evidence) 


Prior 


1 Lens 


-0.24±0.12 


-0.02±0.10 


n CQ7+0.0U2 
"•-'^'-0.038 


17.9±1.6 


476±7 


0.38" 


86/3 


-62 


e <0.6 


Ext. Perturb. 


0.01±0.05 


-0.26±0.08 


0.499±0.04 


2L5±0.9 


447+3 


0.05" 


0.96 


-20 


e <0.6 



Table 1. Parameters of the lens inferred from two optimisations: First line: A single halo models the lens potential. Second line: We add an external 
mass perturbation on top of the halo lens potential. Coordinates are given in arcseconds with respect to the centre of the galaxy deflector, e is the 
central halo ellipticity. Error bars are given for Icr confidence level. 




Fig. 3. Three luminosity maps of the same luminosity Lext at difl'erent smoothing scales as indicated on each panel. We draw in white luminosity 
isodensity contours of 10' and 10^ Lq arcsec"^. We can appreciate how the smoothing scale s influences the shape of the luminosity isodensity 
contours, and by construction the shape of the resulting mass distribution, hence its lensing properties. 



& Kneib (2009), we found that such values of rcut was ensur- 
ing a smooth and extended density profile. We estimate the dPIE 
central velocity dispersions cr, by inverting the equation 

'^z = Mi.j^j (3) 

where is the surface density at the grid nodes location. 
Mij is a mapping matrix whose coefficients depend of the dPIE 
core and cut-off radii (see Julio & Kneib 2009). In order to pre- 
vent negative cP:, we invert equation 3 by a mean-square min- 
imisation technique. Finally, each luminosity map is scaled to 
the corresponding mass map by multiplying all dPIE potential 
velocity dispersions by a constant mass-to-light ratio MgxtA-ext 
. This mass-to-light ratio Mext/Lext is the second free parame- 
ter describing the perturbation produced by the galaxy group. 
The density threshold controls the grid resolution, and as such 
might be considered a critical parameter but we have tried to 
use smaller thresholds, down to 10^ M0 arcsec"^, and the re- 
sults were unchanged. We therefore keep the 10' Mq arcsec"^ 
threshold because the corresponding mass maps require less 
mass clumps. We also force the algorithm to stop after 4 levels of 
splitting. On average, a grid contains about 200 dPIE potentials. 

4.3. Modelling the Lens Accounting for the External 
Perturbation 

We now come back to modelling the SL system, taking into ac- 
count the external mass perturbation parametrised by a smooth- 
ing scale s and a mass-to-light ratio Mext/Lext- We generated 
mass maps with smoothing scales s ranging from 1 to 40" by 
steps of 2.5" and mass-to-light ratios from 10 to 190 by steps 
of 20. Each mass map is then included in the modelling of the 
SL system. This modelling is performed in the image plane. 
Parameters of the lens potential are taken from Section 3.4. For 



each set of parameters (s, Mext/Lext), we quantify the goodness of 
the SL modelling using the image plane RMS, the corresponding 
and the Bayesian Evidence. 

In the following, because our goal is to constrain the galaxy 
group as a whole, we use M/L corresponding to the total mass- 
to-light ratio of the group; i.e. M (resp. L) is the sum of the exter- 
nal mass (luminosity) perturbation and the mass (luminosity) of 
the lens. We checked that degeneracies of each mass component 
near the lens are small. For the range of parameters (s, Mext/Lext) 
investigated in this work, we compute Mext/Miens in a circle of 
radius 10" centred on the lens. This ratio falls between 10 ' and 
10-5. 

Total masses and luminosities are computed within a region 
of 10' X 10' centred on the lens. At the redshift of the group, it 
corresponds to ~ 3 x 3 Mpc^ (Fig. 1). 

5. Results: Properties of SL2S J08544-0121 

For a certain range of parameters characterising the external 
mass perturbation we obtain excellent fits for the observed con- 
straints. We present first the best-fit model, and then the derived 
constraints on the galaxy group properties. We emphasise that 
what we achieve here is to constrain the properties of the galaxy 
group as a whole (on scales of 100") based on a local SL analysis 
only (on scales of 10"). 

5.1. Best-Fit Model for the Lens 

Results are given in Table 1. The best fit is found for a total 
mass-to-light ratio ^ 75 (/ band, solar units, not corrected for 
evolution) and a smoothing scale ^ 20" (Fig. 4). The RMS er- 
ror between observed and modelled image positions in the im- 
age plane is 0.05", yielding a reduced of 0.96. This is a sig- 
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Fig. 4. Results: constraints on the galaxy group mass-to-light ratio M/L 
and smoothing scale s that characterises the size of dark matter clumps. 
Vertical dashed lines are lower-higher constraints (1 cr error bars) on 
group mass-to-light ratio from an independent weak-lensing analysis 
(Section 7). 



nificant improvement compared to the modelling without exter- 
nal mass perturbation, which had RMS =0.38": the difference 
in log(Evidence) values, which accounts fully for the additional 
complexity of the new model (extra parameters and their prior 
PDFs) is 42. The log(Evidence) difference between two mod- 
els compares the relative probabihty of getting a set of data ex- 
plained by the associated model. Here the data set associated 
with the perturbed model is e^^ times more probable than the 
data set associated with the simple model. 

We find the position of the halo to coincide with the centre of 
the light distribution. The modelled position angle of the halo is 
21.5 deg. Comparing this value to the position angle of the light 
distribution is complicated by the bimodal hght distribution of 
the bright galaxy (Section 3.2). In particular, a merger wiU affect 
the light and mass distributions so that agreement may be not 
necessarily expected. 

5.2. Constraints on the Group as a Whole 

differences between each model are translated into confidence 
levels and drawn on Fig. 4. Considering the 2cr contour, we find 
M/L = 98±27 {i band, solar units, not corrected for evolution) 
andi = 20"±9. 

5.3. Is Mass Traced by Light ? 

We are able to reproduce accurately the observational constraints 
when considering an external mass perturbation drawn from the 
light distribution. Because our SL analysis is sensitive to the 
mass, this finding is consistent with the hypothesis that light is a 
good tracer of mass. We note, however, that we have not demon- 
strated the uniqueness of the smoothed light model. 



6. Local Effect of Large Scale Perturbation 

In this Section, we propose to explain why a local SL analysis 
is able to constrain global properties of the galaxy group hosting 
the lens. First, we investigate the lensing impact of the external 



perturbation on the local SL modelUng (i.e. on the local image 
positions). 

6.1. Lensing Properties of the External Perturbation 

The lensing properties of a mass distribution are commonly 
parametrised by a shear y and a convergence k (see, e.g. 
Schneider et al. 1992, for the definition of these quantities). Here 
we estimate the mean shear and convergence experienced locally 
by the lens for each set of parameters (s, Mext/Lext) by averaging 
y and k generated by the perturbation within a 7" square encom- 
passing all the multiple images. 

Fig. 5 shows K and y maps generated by an external per- 
turbation of fixed mass (5.7x101'^ Mo) for different values of s 
(reported on each panel). These maps have been generated for a 
source redshift equal to 1.268. Red crosses indicate the lens po- 
sition. One can appreciate how the experienced shear and con- 
vergence are correlated with the smoothing scale. 

6.2. Shape of the Constraints 

Fig. 6 shows lines of constant k and y on the results of Fig 4. 
We see that the constraints inferred from the SL analysis do not 
follow K lines but do follow y lines. In particular, the best model 
generates a shear of ~ 0.075. We interpret the shape of the con- 
straints as follows: one needs to generate a shear value of ~ 0.075 
with a mass distribution parametrised by a smoothing scale and a 
total mass. For a given total mass, the smoother this mass distri- 
bution (the higher s), the smaller the generated shear. Therefore, 
the smoother the mass distribution, the higher the total mass in 
order to generate a given shear level. 

6.3. Why Closed Contours ? 
Shear, Convergence, and Beyond 

Fig. 6 suggests that the observational constraints require locally 
an external shear component of ~ 0.075. However, our analysis 
rejects some external mass perturbations characterised by such 
a shear (Fig. 4). This suggests, as found by Keeton & Zabludoff 
(2004), that the shear approximation fails to capture all of the 
environmental effects. In other words, the shape of the con- 
straints follow constant shear lines. However, the contours do 
close, which means that the constraints are sensitive to more than 
the shear; most probably higher order terms beyond shear that 
are naturally provided by the modelUng proposed in this work. 

Besides, we estimate the shear experienced by images 1.1, 
1.2, 1.3 and 1.4. To do so, we consider all models falling in the 
1<T contour. For each model, we compute the shear experienced 
by the images, and from these numbers, we estimate the mean 
shear and the associated standard deviation. They are equal to 
0.075, 0.074, 0.073 and 0.073 respectively, with a typical un- 
certainty equal to 0.01. Therefore, each image does experience 
the same shear value within the error bars. We will investigate 
further the differences between our approach and a constant ex- 
ternal shear approach further in Appendix. 

7. Mass-to-Light Ratio from Weak Lensing 

We have presented constraints on the mass-to-light ratio of 
galaxy group SL2S J08544-0121 based on a local SL analy- 
sis. In this Section, we aim to constrain the mass-to-light ratio 
of SL2S J08544-0121 from an independent weak lensing (WL) 
analysis, which is intrinsically more sensitive to the projected 
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Fig. 5. Shear (upper panels) and convergence (lower panels) maps generated by the external perturbation and experienced by the SL system whose 
centre is given by the red crosses. The total mass is fixed to the same value (5.7 lO'"* Mq) in all panels. Panel sizes are 600 x 600 square arcseconds, 
and the smoothing scales j vary as indicated on each panel. White contours correspond to shear levels of 0. 1 and convergence levels of 0. 1 and 
0.2. One can appreciate how the shear and convergence generated by the group are correlated with the smoothing scale. 



mass distribution on large scales. The goal is to check whether 
M/L inferred from the two techniques are consistent. 

For a detailed description of our WL methodology, see 
Limousin et al. (2009). Here we just give a brief reminder. We 
select as background sources all galaxies with /-band magnitudes 
in the range 21.5-24. Their density is 13.5 arcmin"-. The com- 
pleteness magnitude in this band is 23.91. The seeing is 0.51". A 
Bayesian method, implemented in the Im2shape software (Bridle 
et al. 2002), is used to fit the shape parameters of the faint back- 
ground galaxies and to correct for PSF smearing. From the cata- 
logue of background galaxies, Limousin et al. (2009) performed 
a one-dimensional WL analysis. They fit a Singular Isothermal 
Sphere model (SIS) to the reduced shear signal between 150 kpc 
and 1.2 Mpc from the group centre, finding an Einstein radius of 
5.4±2.1". In order to relate the strength of the WL signal to a 
physical velocity dispersion characterising the group potential, 
Limousin et al. (2009) estimate the mean geometrical factor us- 
ing the photometric redshift catalogue from the T0004 release 
of the CFHTLS-Deep survey ^ (lenna & Pello 2006). They find 
csis = 658;^^^ kms"'. This translates into a total mass within 
the considered square of 5.3±2.0 lO'** Mq. Because the total lu- 
minosity is 5 10'^ L0, we find a mass-to-light ratio equals to 
106±40 (i band, solar units, not corrected for passive evolution). 



http://www.ast.obs-mip.fr/users/roser/CFHTLS_T0004/ 



This is comparable to the M/L constrained by SL only. The 
good agreement (Fig. 4) between the two methods gives support 
to the SL only analysis. 



8. Discussion 

8.1. Mass is Traced by Light 

An external mass perturbation derived from the light of the group 
members allows us to fit accurately the observed constraints. 
Because observed constraints are sensitive to mass and not to 
light, this suggests that light is a good tracer of mass. We note 
that this assumption constitutes an efficient way of taking into 
account an external mass perturbation in SL modelling. Indeed, 
this perturbation is fully described with only two parameters, the 
mass-to-light ratio and the smoothing scale. In contrast, describ- 
ing this perturbation parametrically using a mass clump would 
require at least three parameters (position and velocity disper- 
sion), unless independent data motivate additional priors (see, 
e.g. Tu et al. 2009, where X-ray observations allow one to con- 
strain the group centre). 

8.2. Wliat is fhe source of tine Constraints ? 

Why our SL analysis did allow us to infer some constraints on 
the whole galaxy group? We claim that this is due to the per- 
turbed state of the SL system of SL2S J08544-0121. Most of the 
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Fig. 6. As in Fig. 4, constraints obtained on the galaxy group as a whole derived from the local SL analysis are shown as dashed contours. Solid 
lines corresponds to lines of constant k (left) and y (right) generated by the external mass perturbation and experienced locally by the lens, then- 
values are labelled on each line. 



perttirbed signal of the multiply imaged system comes from im- 
age 1 .4, because it is located much further from the lens centre 
(~ 8") than images 1.1, 1.2 and 1.3 (~ 5"). If we remove image 
1.4 from the set of observational constraints, we are able to fit 
very well the remaining images without considering any external 
mass perturbation (the lens being modelled as in Section 3.4.) In 
that case, we get RMS = 0.03" and a reduced equals to 0.03. 
Therefore, ignoring image 1.4 prevents us to put any constraints 
on the external mass perturbation, i.e. the host galaxy group. 
This shows that image 1.4 yields the constraints presented in 
this work. This finding will help us diagnose the type of the SL 
systems favouring this kind of analysis (see Section 8.5). 

We note that the SL analysis presented here is very simple 
since we just conjugate a couple of images with each other. In 
particular, we do not use the constraints coming from the whole 
Einstein ring. More sophisticated methods fuUy take into ac- 
count arc brightness constraints (see, e.g. Warren & Dye 2003; 
Suyu et al. 2006; Barnabe & Koopmans 2007). We are aware that 
we ignore some informations that could allow us to put stronger 
constraints on the galaxy group. One the other hand, the basic 
level of the SL analysis done here emphasises even more the 
prospects of this method. 

8.3. The Satellite Galaxy 

We have assumed that the satellite galaxy does not produce sig- 
nificant shear on the images used as constraints. However, one 
could argue that neglecting this satelhte galaxy effectively pro- 
duces the claimed constraints from the SL analysis. This is not 
likely, because of the location of the satellite galaxy with respect 
to the multiple images (Fig. 1). The dwarf galaxy may produce a 
marginal shear on on images 1.1 1.2 1.3, but is unlikely to have 
any significant influence on image 1.4, the one image yielding 
most of the constraints. Indeed, the distance between the satellite 
galaxy and image 1.4 is ~ 13". We note that we do not quantify 
the bias that could result under our working assumption. 

Besides, a paper focusing on the properties of the satellite 
galaxy (Suyu & Halkola, submitted) shows that even with the 
satellite galaxy included in the lens model, an external shear of 



approximately the same magnitude is needed to fit the observed 
constraints. 

8.4. Choice of the Lens' Scale Radius 

The dPIE scale radius is where the logarithmic slope of the 3D 
density profile smoothly decreases from -2 to -4. The scale ra- 
dius of the lens is set to 250 kpc in the present analysis. We have 
also done a complete analysis for a scale radius of 400 kpc as a 
sanity check and found that the results inferred for the group do 
not change significantly. To understand why, we superimposed 
critical lines of the best- fit parameters of Table 1, for a source 
redshift of 1 .268 (without extemal perttirbation), and the critical 
lines of the best-fit parameters of Section 5.1 (with external mass 
perturbation). We find that the external mass perturbation gener- 
ates a critical line shift of 1.3". In parallel, we investigated the 
critical Unes shifts between various scale radii; increasing from 
250 kpc to 400 kpc and decreasing from 250 kpc to 100 kpc. The 
shifts are 0. 12", an order of magnitude smaller than the shift due 
to the extemal mass perturbation. 

8.5. Looking for Perturbed SL Systems 

We propose to characterise the kind of perturbed SL systems one 
should look for to be able to perform analyses similar to the one 
presented in this work. From the ring test done in Abell 1689 
(Tu et al. 2008, see Section 1.3) and the analysis presented in 
this paper, we hint the need for a non-dominant SL system used 
as a particle test probing the main potential. 

This is linked to the global geometry of the structure hosting 
the SL system: in order to be perturbed, a SL system should not 
be at the centre of the structure. Indeed, if the lens studied in this 
paper would have dominated the whole group potential, image 
1 .4 would have been located at a similar distance from the lens 
centre as images 1.1, 1.2 and 1.3. 

The Cosmic Horseshoe (Belokurov et al. 2007; Dye et al. 
2008) illustrates this point: it is an almost complete Einstein ring 
of radius 5" containing a luminous red galaxy in its centre. This 
galaxy is the brightest object in the group of ~ 26 members as 
revealed by the SDSS photometry and it dominates the group 
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light distribution. No external shear is required in the model of 
the Cosmic Horseshoe SL system, which is already suggested by 
the nearly perfect circle outlined by the ring. 

To summarise, we should look for multiply imaged systems 
where one of the images is found at a larger radial distance than 
the other images of the SL system. 

8.6. SL2SJ08544-0121: Further Evidences for a Bi modal 
Mass Distribution from Spectroscopy of Group Members 

We have shown that the strong lensing modelling only provides 
strong hints for a bimodal mass distribution: a first one is clearly 
associated with the strong lensing deflector, and the second one 
perturbing the strong lensing configuration seems, to first order, 
associated with the second light peak of the bimodal Ught map. 
This suggests a dynamically young structure in the process of 
formation. Additional data support this hint. We measured red- 
shift for 36 galaxies along the direction of SL2S J08544-0121 by 
using spectroscopic data acquired with FORS2 at the ESO Very 
Large Telescope (VLT), and confirmed the presence of a high 
concentration of galaxies at z ~ 0.35 (Munoz et al., in prepara- 
tion). A careful analysis of the redshift distribution of galaxies 
around this peak reveals two close structures with a radial veloc- 
ity difference of Vr = 1 180kms"\ maybe boimded-incoming 
structures. This result is in agreement with our strong lensing 
only hints. 

9. Conclusion 

We propose a method to constrain the dark matter distribution 
of galaxy groups and clusters. Exploiting information contained 
in perturbed SL systems, we use the SL geometry to probe the 
main potential of the host structure responsible of that perturbed 
state. 

We show that the SL only constraints on the mass-to-Ught 
ratio of SL2S J08544-0121 are in good agreement with WL 
constraints obtained independently, comforting the reliability of 
the proposed method. Moreover, the SL only analysis provides 
strong hints for a bimodal mass distribution, which is confirmed 
by spectroscopy of galaxy group members. 

We advocate the need for a dedicated search of perturbed SL 
systems in the HST archive in order to test and validate fiu'ther 
this method, particularly promising in the light of future large 
surveys that will yield thousands of SL systems, some of them 
being perturbed enough to perform the test presented in this pa- 
per. 
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Appendix A: Taking an External Mass Perturbation 
Into Account: Comparison with Other 
Approaches 

We have proposed in this article a way of taking into account an 
external mass perturbation in a strong lensing (SL) modelling. 
Here we try other possible and more conventional approaches: 
i) a constant external shear profile and ii) a Singular Isothermal 
Sphere centred on the second high luminosity peak, which, by 
construction, is the main mass concentration perturbing the SL 
in the method proposed in this work. 

A.1. A Constant External Shear 

Although unphysical (any mass distribution will not generate 
shear only but also convergence), the external shear model 
is widely used and is often a good approximation. Here we 
address the modelling of the SL system with a constant ex- 
ternal shear component parametrised by a position angle and 
a strength (yKst)- This modelling is performed in the image 
plane. Parameters of the potential describing the lens are set 
as in Section 3.4. The external shear strength is allowed to 
vary between and 0.3. The upper limit corresponds to a very 
strong shear value: for comparison, the massive galaxy cluster 
Abell 1689 (Limousin et al. 2007), produces an average shear 
value of 0.23 at 50" away from its centre. 

We are able to get a very good fit, with < ^- The best 
model corresponds to a circular halo for the lens (e = 0.036) 
centred on the bright galaxy (X - -0.36", Y = 0.07"), making its 
position angle (95 degrees) irrelevant. The lens fiducial veloc- 
ity dispersion equals to 450+8 kms ' (lo"). The external shear 
is described by yKst = 0.19, more than twice the one derived in 
Section 6.2 (~ 0.075), and a position angle equal to 19.8 degrees. 

The Lenstool software does explore the parameter space us- 
ing a MCMC sampler (Julio et al. 2007). Therefore, we can 
use these MCMC realisations in order to investigate the degen- 
eracies between the different parameters. The following figures 
have been generated this way. 

Fig. A.l shows that there is a strong degeneracy between e 
and yKst- We see that the solution derived in Section 5.1 (i.e. 
an external shear of ~ 0.075 and an ellipticity of ~0.5) is in- 
cluded in the 1 -cr contour. On the other hand, the position angle 
of the external shear is very well constrained to be ~ 20 degrees. 
This position angle points towards the second high luminosity 
light clump. This suggests that, to first order, the external mass 
perturbation is dominated by this component. We note that the 
best-fit model needs an external shear of order 0.18, which is a 
pretty unlikely value in our case because it is comparable to what 
would be experienced at ~ 100" from the centre of Abell 1689. 

A.2. An SIS profile 

The first order mass perturbation is associated to a second peak 
of high luminosity. We put an SIS mass distribution at the lo- 
cation of this second luminosity peak (X,Y = -53,10" wrt the 
lens). We allow its velocity dispersion to vary up to 800 kms"', 
an upper limit motivated by the WL analysis of the full group 
(Section 7), and do the SL modelling with parameters set as in 
Section 3.4. We are able to get a very good fit, with;^^^ < 1. The 
lens halo is centred on the bright galaxy. Its ellipticity equals 
0.43^g p and its position angle 27+2 deg. The lens fiducial ve- 
locity dispersion equals to 441+7 kms"' (lo"). The external 
shear and convergence generated by the SIS profile at the lo- 
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Fig. A.l. Results of a constant external shear model. Top: degeneracy 
between the halo ellipticity e and the strength of the external shear yKst- 
Bottom: posterior probability distribution for the position angle of the 
external shear 



cation of the multiple images are equal by definition, between 
0.07 and 0.10. 

We show on Fig. A.2 the degeneracies between the lens halo 
ellipticity and the SIS profile velocity dispersion, related to the 
strength of the external shear experienced by the multiple im- 
ages. We see that the solution derived in Section 5.1 (i.e. an ex- 
ternal shear of ~0.075, corresponding to osis ~ 700 kms"' and 
an ellipticity of ~0.5) is included in the 1-cr contour. 

A.3. Discussion 

In each cases investigated in this Appendix, we find that the solu- 
tion we have derived in Section 5.1 using our original method is 
consistent with solutions derived with more conventional meth- 
ods. 

We note that conventional methods exhibit strong degenera- 
cies between the lens halo ellipticity and the strength of the ex- 
ternal shear. These degeneracies are smaller in the case of the 
SIS profile (lens ellipticity is constrained between 0.3 and 0.6) 
compared to the case of a constant shear profile (ellipticity un- 
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Fig. A.2. Degeneracies between tiie lens tialo ellipticity and tlie SIS pro- 
file velocity dispersion, related to the strength of the external shear ex- 
perienced by the multiple images. 



constrained between the allowed priors: and 0.6). The main 
difference between the SIS profile and the constant shear profile 
is that the SIS profile generates both shear and convergence. 

Concerning the lens itself, we note that all fitted fiducial ve- 
locity dispersions are consistent, whatever the method used to 
take into account the external mass perturbation. They fall be- 
tween 433 kms"' and 458 kms"'. This translates into a pro- 
jected mass computed in a radius of 10" between 0.93 and 
1.04x10'-' M0. This is expected because the mass of the lens 
within this radius is set by the location of the SL constraints and 
therefore does not depend much on the external mass perturba- 
tion (see also Section 4.3). 



